Reentrance of Bose-Einstein condensation in spinor atomic gases in magnetic field 
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We calculate the Bose-Einstein condensation (BEC) temperature of spin-1 atomic bosons in ex- 
ternal magnetic field, taking into account the influence of the quadratic Zeeman effect. In case that 
the quadratic Zeeman coefficient is positive, the BEC temperature exhibits a nontrivial dependance 
on the magnetic field and a magnetic-field-induced reentrant phenomenon of BEC is observed. This 
phenomenon could be well understood by the competition between the linear and quadratic Zeeman 
effects. Reentrance of BEC in a trapped spinor Bose gas is also discussed. 
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I. INTRODUCTION 

Since the first observation of spin-1 ■^■^Na Bose-Einstein 
condensation (BEC) in the optical trap [110], Bose gas 
with spin degrees of freedom, called the spinor Bose gas 
d, 3 , becomes one of the central topics in cold atomic 
physics. The spin degrees of freedom permit a variety of 
additional exotic phenomena to be manifest. Consider- 
able experimental and theoretical works have been ded- 
icated to investigating static or dynamic properties rel- 
evant to spins of the condensate, including magnetic or- 
dering [1,13,00, coherent spi n-dy namics 
spin-domains or textures [l2, Il3l . [3, [Tsl . 
on. In addition, a number of researchers have laid eyes on 
finite-temperature properties of spinor bosons, especially 
on the BEC phenomenon itself [ll, [H, HO, [HI, HI [H . 

Bose-Einstein condensation in a free spin-1 Bose gas 
has ever been studied by Yamada [ll|, Simkin and Co- 
hen [l^. The obtained results suggest that the Bose 
gas is more sensitive to external magnetic field than a 
Fermi gas and the BEC temperature rises up with the 
magnetic field due to the Zeeman effect. If there are 
ferromagnetic (FM) couplings between atoms, BEC oc- 
curs in a more fascinating way because the FM coupling 
induces a new phase transition, the FM transition ^2Q] . 
Its critical temperature is called the Curie point. It is 
already demonstrated that the Curie point is never be- 
low the BEC temperature and both critical points are 
increased by the FM couphng [g^]. These conclusions 
have also been confirmed by other groups [2ll|. More- 
over, Bose-Einstein condensation in an optically trapped 
spinor Bose gas has been elaborately studied [22, [2^ . 

In the present study, we concentrate on the Bose- 
Einstein condensation in an = 1 spinor Bose gas by 
considering quadratic Zeeman effect (QZE). As Stenger 
et al. indicated [3|, the Zeeman energy of a spin-1 Bose 
atom is given by 



El^ = Ea- pha + qh^a^ 



(1) 



where a refers to the spin-2; index of hyperfine state 



\F = 1, mp = a) [a ^ +1, 0, —1), h to the magnetic field, 
and Eq to the Zeeman energy of the mp — state. The 
second term represents the linear Zeeman splitting and 
the last term arises from the quadratic Zeeman effect. 
For a spinor atom, the quadratic Zeeman effect can not 
be neglected due to the hyperfine structure [13, |2^ . Very 
recently it is found that the QZE can even be induced by 
the dipole trap, which is used to confine cold atoms ex- 
perimentally, for the spinor atom ^^Cr [25j . 

The QZE plays an important role for the understand- 
ing of spinor atomic bosons, especially in strong magnetic 
field. Various of effects caused by the quadratic Zeeman 
shift have been reported in recent years. It, together 
with the linear Zeeman effect and spin-dependant inter- 
actions, can not only lead to a number of novel ground 
states [3, [2^, but also affects significantly on spin dynam- 
ics [3, [J of the spinor condensates. The spin dynamics is 
dramatically suppressed in magnetic field, owing to the 
quadratic Zeeman splitting. The QZE can also influence 
on the vortex state, and bring about new type vortex 
states in spinor condensates j26|. In this paper, we aim 
at discussing the Bose-Einstein condensation of the spin- 
1 Bose gas. Interestingly, we show that the QZE can 
induce reentrance of BEC in the Bose gas. 

The spin-1 Bose gas is described by the following 
Hamiltonian: 



n 



dr 



-phi^iF:,^j,+qh^iji{F:,fi^t, 



(2) 



'Electronic address: ' qgu@sas.ustb.edu.cn| 



where ipa(j) is the field annihilation operator for an atom 
in hyperfine state |l,a) at point r. The first term rep- 
resents kinetic energy. The terms with coefficients cq 
and C2 denote the spin-independent and spin-dependent 
interactions, respectively. The forth and fifth terms de- 
scribe the linear and quadratic Zeeman effects, respec- 
tively. We suppose that the linear Zeeman coefficient p 
is always positive hereinafter, but the quadratic Zeeman 
coefficient q can be either positive or negative. F is a 
vector which consists of three components of the 3x3 
Pauli spin-1 matrices [1, [3|- 
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This paper is organized as follows. In Sec. II, we cal- 
culate the condensation temperature and the condensate 
fraction of homogenous spinor gases in the external mag- 
netic field. The reentrance phenomenon is observed in 
case of (J > and a detail explanation is given on how 
it is induced by the quadratic Zeeman effect in a free 
gas. The Sec. Ill extends the study to the trapped inter- 
acting gases. The interaction between atoms is treated 
by using a mean-field approach. It is suggested that the 
interaction tends to decrease the condensation tempera- 
ture, but the reentrance phenomenon sustains. The last 
section gives a brief summary. 
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II. REENTRANCE IN A HOMOGENOUS GAS 

First we consider a free spin-1 Bose gas, ignoring the 
interactions in the Hamiltonian Then the effective 
Hamiltonian for the grand canonical ensemble reads 



H - Nfi = ^ (£fe - ^ - pha + qh^a"^ 



)nko 



(3) 



where — fi?k^ /2m* is the kinetic energy of free parti- 
cles with mass to*, ^ is the chemical potential, and N is 
the operator of total particle number. Since the Hamil- 
tonian is diagonal, we may calculate the grand thermo- 
dynamical potential 



n = --\nZ = --InTr 



(4) 



where Z = Tr 



is the partition function, and 



(3 = l/ksT. The density of particles is derived as 



1 / an 

V \d^J.a/ Ty 



where V is the volume of the system, and n = N/V is the 
particle density. From Eq. ([5]), we can obtain the basic 
equation determining the phase diagram of the spin-1 
Bose system, 



1 



+ /3/2 



m*kBT 
27rfi2n2/3 



3/2 



/3/2 



^ + ph — qh? 



+ /3/2 



keT 
fi — ph — qh^ 



(6) 



where ric = fic/n is the condensate fraction, he is the con- 
densate density, and is the polylogarithm function 
defined as [201 



oo / x\k 



(7) 



k=l 



Above the condensation temperature T^, the conden- 
sate fraction Uc — and the chemical potential fi is de- 
termined by Eq. ([5]) for each given temperature. At Tc, 



FIG. 1: Reduced BEC temperature tc vs external magnetic 
field h for the spin-1 Bose gas with g = 1. The inset plots 
BEC temperatures for all the three cases of g = — 1,0 and 1 
so as to make a comparison. 

atoms lying on the lowest energy sub-level begin to con- 
dense. In order to determine Tc, one have to make it 
clear which sub- level is the lowest first. From Eq. ^ it 
is easy to see that the energy level for the itif = I state is 
lower than the other two as long as the external magnetic 
field h is sufficiently weak, so the nip = ^ atoms condense 
once the temperature is below T^. The chemical potential 
fi satisfies the equation ^ + ph — qh? < above Tc, and 
fi + ph — qh'^ = at and below Tc- Then the condensate 
fraction is calculated by the following equation 



/3/2 



T 

c 



3/2 



1 



3C(3/2) 
qh^ 



C(3/2) 



(8) 



where 



T° — 



m*kB V3C(3/2) 



2/3 



is the BEC temperature for a free spin-1 Bose gas with- 
out external field and C(3/2) = f 3/2(0) w 2.612 is the 
Riemann zeta function. 

As the external magnetic field becomes stronger, the 
energy level structure may be modified because of the 
influence of the QZE, especially in the case of g > 0. 
When qh^ > ph, the mp — sub-level becomes lower 
than the mp — ^ sub- level, and tof = atoms condense 
in place of mp = 1 atoms correspondingly. Therefore, 

^ as T ^ Tc. Below Tc, Eq. © reduces to 



1 



T 



3/2 



1 



-C(3/2) + /3/2 



3C(3/2) 
/ ph 



/3/2 



ph ~ qh^ 
knT 



V keT 



(9) 



Based on above discussions, we produce the tc-h phase 
diagram numerically, as plotted in Fig. [TJ In our calcu- 
lations, a reduced temperature t = T/[r^"(3C(3/2))2/3] = 



3 




a) 



FIG. 2: Condensate fraction ric vs external magnetic field h 
at different temperatures t = 0.2, 0.3, 0.4. 



{m*kBT)/{2'Kh^'n?^^) is introduced, and thenph and qh? 
must be re-scaled by the factor of fcBr^"(3C(3/2))2/3. The 
particle density is suppose to be 1; p, q and h are cho- 
sen to be dimensionless parameters and by setting the 
re-scale factor to be 1. For simplicity, we fix the value of 
p as p = 1 and place emphasis on discussing three dis- 
tinct cases with respect to the QZE, q = —1, 0, and 1. As 
shown in the inset of Fig. [1] the phase diagram for the 
g = — 1 or case is very simple. The BEC temperature 
tc grows monotonously as h increases. In contrast, the 
q = 1 case is very intriguing: tc increase with h first, and 
then decreases. After dropping to a minimum value at 
h — 1.0, it grows up again. During the temperature range 
oi t ^ 0.33 io t 0.43, a typical reentrant phenomenon 
of BEC induced by the magnetic field h is observed. 

To further demonstrate the reentrant phenomenon in 
spin-1 atomic bosons with g = 1, we calculate the con- 
densate fraction ric at some given temperatures, as shown 
in Fig. [21 The condensed particles become dramatically 
less around h = 1.0, corresponding to the minimum point 
of tc in Fig. [1] Especially at i = 0.4, the condensate van- 
ishes around h = 1.0, and then reappears in stronger 
field. The reentrant phenomenon can be interpreted as 
the consequence of the competition between linear and 
quadratic Zeeman effects, as discussed in the following. 

Without magnetic field, the energy spectra for the 
three hyperfine spin states, = Ek—fJ.—pha+qh'^a^, are 
degenerate. Therefore, each spin state is equally occu- 
pied, fii = no = n_i = n/3. n/3 is just the critical par- 
ticle density, which determines the BEC temperature to 
be <2 = {l/[3C(3/2)]}2/3 K 0.2535. When external mag- 
netic field is applied, hyperfine energy levels are shifted 
due to both the linear and quadratic Zeeman effects. Fig. 
13] shows schematically the splitting of energy levels for 
atoms with positive quadratic Zeeman coefficient, q > 0. 
The energy dispersion with respect to k is not depicted 
explicitly. 

When h is relatively week, the linear Zeeman effect 
dominates over the QZE and the energy level splitting 
is shown in Fig. ^a), which satisfies E\ < E'^ < E^^. 
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£i^- H + qh^ - ph 



Sf.- /! + qh^ + ph 
£i^- H + qh^ 
Si^— fi + qh^ - ph 



FIG. 3: Energy splittings of the three spin components in case 
of g > with h being relatively weak (a) or strong (b). 



Consequently, ni > ng > f^-i, and fii acts as the crit- 
ical density of particles. Because fii > n/3, the BEC 
temperature is larger than i^. As h increases gradually, 
the three energy levels become more separated so that ni 
and thus tc tend to growing. This region of h is labeled 
as I in Fig. [1] When h = p/{2q) = 1/2, the separation 
between El and E^ arrives at the largest value, p^ /{4,q). 
However, Fig. [l]shows that Region I terminates at a field 
h\ larger than 1/2. During 1/2 < h < h\, though the 
separation between E], and E^ becomes shrinking, the 
separation between E\ {E^) and E^^ is still enlarged. 
As a result, more atoms occupy the = and mp = 1 
so that fii keeps increasing with /i, until h = h*. 

In Region II of Fig. [U the quadratic Zeeman term qh^ 
becomes dominating, with the result that E^ becomes 
closer to the E^ significantly. Although E^ is still the 
lowest energy level, the occupation number ni begins to 
decrease and tc drops accordingly. But otherwise, no in- 
creases. Since E^^^ is far higher than E'^ and El, few 
atoms occupy the mp — —I energy level. When h ap- 
proaches a particular value, /ij — p/q = 1, E^ = E\ and 
hence no = ni which is slightly smaller than n/2. So 
the BEC temperature tc ~ 0.3327 is slightly higher than 
{■i/2f/Hl w 0.3321. 

Increasing h further, E'^ becomes the lowest energy 
level in place of E\, as indicated in Fig. [3jb). Therefore 
the condensed atoms are in the mp — Q state, instead of 
in the nij;- = 1 state. Meanwhile the BEC temperature 
tc grows up again. As the stronger h makes the three 
energy levels separate farther and farther, no — > n and 
tc i^^^tf) in the high field limit. This point is noted in 
Region III of Fig. [T] 

In case of g = or g 
est energy level all the way in magnetic field and only 
the nip — 1 atoms can condense. The magnetic field 
h persists in enhancing splittings of the three energy 
levels. Correspondingly, the BEC temperature grows 
monotonously with h tending to stronger, as shown in 
the inset of Fig. [TJ 



'1, El. remains the low- 
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FIG. 4: Reduced BEG temperature tc vs external magnetic 
field h for a ferromagnetically coupled spin-1 Bose gas. The 
coupling |c2| = 0,0.2,0.4,0.6,0.8, and 1 from bottom to top. 



We note that the reentrant phenomenon has been dis- 
covered in a number of quantum many-particle systems, 
especially in some unconventional superconductors. For 
example, it is reported that there appear successively 
two distinct superconducting phases with the pressure 
in heavy fermion material CeCu2Si2 [27'|. A magnetic- 
field induced reentrance of superconductivity has ever 
been observed in a pseudo-ternary Eu-Sn material [28j . 
For cold atoms, Kleinert et al. have discussed reentrance 
in the quantum phase transitions of an interacting Bose 
gas d^. They suggested that Tc shifts upward with the 
interacting strength first, then it is suppressed to zero. 

The *^Rb gas happens to exhibit a positive QZE 
and therefore it may be an appropriate example where 
the reentrant phenomenon is expected to be observed. 
There exist weak FM couplings between ®^Rb atoms, as 
described by the C2 term in Eq. ^ with C2 < 0. So 
we proceed to examine how the FM coupling affects on 
the reentrance of BEC, by decoupling the C2 term via 
me an- field theory |20|, 



= -|c2|m^(nfci -hk-i), 



(10) 



where m is the normalized magnetization defined as m = 
X]fe("'fci ~ nk-ij/N. Below the BEC temperature, it is 
given by 



m 



3/2 



1 



T^J 3C(3/2) 



C(3/2) - /3/2 



2ph + 2|c2|m 



(11) 



if condensed atoms are in the mp — \ state or 



T 

-/s/; 



3/2 



3C(3/2) 



/3/2 



ph + |c2|m — qh? 



knT 



f ph+ \c2\m + qK^ 



if condensed atoms in the mp — ^ state. 



(12) 



Based on Eq. (|TC)|) . when the C2 term is included, Eqs. 
© , dH]) and © should be rewritten by replacing ph with 
ph-\-\c2\m. Hence it seems that the FM coupling performs 
a role to enhance the linear Zeeman effect. Combining 
Eqs. dH) and flT]), or ^ and ([12]), we obtain that the 
reentrant phenomenon still manifests, as seen in Fig. 4. 
The FM coupling tends to stabilize the mp = 1 conden- 
sate, so that the BEC temperature is upraised in the weak 
field region oi h < with the FM coupling increasing. 
Meanwhile, moves to the high-field side. 



III. REENTRANCE IN A TRAPPED 
INTERACTING BOSE GAS 

In experiments, atomic gases are usually confined in an 
effective three-dimensional harmonic trap, so it is useful 
to extend the above calculations to the trapped case. In 
this section, we consider a system of interacting atoms 
trapped by the following potential 



V{v) 



1 



where uji(i = x, y, z) denotes the angular frequency. 

There are two interaction terms in the Hamiltonian 
Actually, the spin-dependent interaction is rather weak 
in comparison with the spin-independent interaction in 
cold atoms. For example, in the ^^Rb gas C2 is about 2 
orders of magnitude smaller than cq. Therefore, in the 
following calculations we neglect the C2 term and turn to 
deal with the cq term in some detail. 

For a trapped scalar Bose gas, the interaction can be 
well-treated by the mean-field theory. The obtained re- 
sults on the BEC temperature and the condensate frac- 
tion are quite accurate [13, HH . The BEC temperature is 
shifted to lower temperatures in the presence of repulsive 
interactions [30]. Hereinafter, we calculate the BEC tem- 
perature and condensate fraction for spinor Bose gases in 
the magnetic field based on a generalized mean-field ap- 
proach from the scalar to the spinor Bose gases [23j . 

According to the mean-field theory [3l[, the field op- 
erator could be treated as tpa = 4'<t + Sipa- with 0o- = 
('0cr). Within the Hartree-Fock (HF) approximation, the 
condensate can be described by the finite-temperature 
generalization of the spinor Gross-Pitaevskii equations 
(GPEs) nil, which are given by 



ih 



d 



2m* 



+ V{r) 



— pha + qh'^ + CQ{n + n^)j (jja- (13) 
for the mp = a = ±1 condensation and 

o ?-2 

^ ( - + W + co(^ + ^0 )) ^0 (14) 

for the mp — condensation. Here nj^ = It/io-p, ri^ — 
{6ip^^6ipa) and n^. — n'^+n'^ describe the condensed, non- 
condensed and total density of the mp = a component. 
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FIG. 5: Reduced BEG temperature tc vs external magnetic 
field h for difi'erent interaction strength. The rescaled inter- 
action parameter = 0,0.05, and 0.1 respectively. 



respectively, n = Ua denotes the total density of all 
atoms. As discussed in Sec. II, the mp = —'^ atoms will 
never condense in the cases we considered. So for Eq. 
(fT3)) . we only discuss a — 1 case. 

The dynamics of noncondensed atoms satisfy the fol- 
lowing equations, 

-V2 + F(r) 



at 



2m* 



ph + qh^ + co(n + ni)jSipi, 

-— V2 + T/(r)+co(n + 7io)UV'o, 

V 2m* / 

V 2m* ^ ' 

+ ph + qh^ +co{n + n^i)jSijj^i. (15) 

Then we derive the effective energy spectra for all the 
three components. 



ot 



ei(k,r) = 
£o(k,r) = 
£_i(k,r) = 



■21,2 



2m* 
2m* 
2m* 



+ V{r) -ph + qh^ + co(n -f rii), 
+ y(r) + co(n + no), (16) 
+ V{v) +ph + qh^ + co(n + n_i). 



Using the semiclassical approximation, the density of 
thermal atoms with mp = a becomes 



e^(/c,r) 



, (17) 



where the chemical potential /i could be derived from 
Eq. p3)) or p4)l . On the basis of the Thomas- Fermi 
approximation, one gets that 

^^V{Y)-ph + qh'^ +co{n + nl) (18) 

when the linear Zeeman effect dominates and only the 
mp = (7=1 atoms can condense, or 
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FIG. 6: (a) Condensate fraction of noninteracting Bose gas in 
harmonica trap at different temperatures t = 0.5, 0.6, 0.7. (b) 
Gondensate fraction of interacting Bose gas in harmonic trap 
at different temperatures t = 0.4, 0.5, 0.6, for Cq =0.1 case. 



at sufficiently large field when the quadratic Zeeman ef- 
fect dominates and the mp — Q atoms become condensed 
instead. 

For simplicity, in our calculations all the energy pa- 
rameters, such as /cbT and cpiV, are rescaled by a factor 
keT* = h\^/^uj^ujyUj,N\^l^ = [3C(3)yV2]'/'fcBT0 
where TO = fi(w^W2,w^iV)i/3/[fcs(3C(3))i/3] is the BEG 
temperature for N trapped free spin-1 atoms without 
external magnetic field. Then the reduced tempera- 
ture reads t = T/T* = \X{i)^JVl2\-^I^T jT"^ and = 

[3C(3)yV2]"^/^ « 0.6048 for the cq = /i = case. 

The BEG temperature can be obtained self- 
consistently from Eqs. and ([TT]), as shown in Fig. 

[5] Obviously, the reentrance phenomenon sustains in 
the presence of the trapping potentials. The cq term 
results in decrease of the BEG temperature at all mag- 
netic fields. As studied previously, interactions tend to 
suppress the BEG in scalar bosons [30| and the BEG tem- 
perature drops linearly with the interaction strength cg 
when Co is small enough. We estimate that the linear re- 
lationship also holds in our case when both the resacled 
interaction parameter Cg and the resacled magnetic field 
h are smaller than 0.01. But in Fig. [5l Cq is up to 0.1 
and thus the linear relationship becomes invalid. 

At last, we calculate the condensate densities below the 
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critical temperatures according to Eq. ((T3)) (at low mag- 
netic fields) or Eq. (at high magnetic region). The 
densities of noncondensed atoms can still been obtained 
from Eqs. and ^T7\ . We get that the condensate 

mainly exists in the center of the trap, while the non- 
condensed atoms mainly appear on the edge region. The 
condensate fractions for both the Cq = and 0.1 cases are 
plotted in Fig. [S]^a) and (b) respectively. Again the reen- 
trance phenomenon is illustrated for both cases. There 
is no much more substantial changes in comparison with 
the homogeneous case. 

IV. SUMMARY 

In summary, we have investigated the influence of the 
QZE on the BEC temperature. In the case of free spinor 
Bose gas which exhibits a positive QZE {q > 0), the BEC 



temperature Tc increases with the external magnetic field 
h first, and then begins decreasing until dropping to a 
minimum value. Tc increases again in stronger h. Thus 
the system shows a reentrant phenomenon of BEC with 
respect h in the phase diagram. Similar properties are 
also expected in trapped interacting Bose gases. We sug- 
gest that this phenomenon can be well understood by the 
competition between the linear and quadratic Zeeman ef- 
fects. We have also calculated the BEC temperature in 
the case of negative QZE {q < 0) and in the presence of 
interactions. 
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